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Abstract 

In this note we carry out the counting of states for a black hole in 
loop quantum gravity, however assuming an equidistant area spec- 
trum. We find that this toy-model is exactly solvable, and we show 
that its behavior is very similar to that of the correct model. Thus 
this toy-model can be used as a nice and simplifying 'laboratory' for 
questions about the full theory. 



1 Introduction 

The present paper is concerned with the description of black holes and cal- 
culation of their entropy in the framework of loop quantum gravity (see 
IH [2) for some general introduction to loop quantum gravity). The liter- 
ature on this subject is large and includes (but is by no means limited to) 
the pioneering work [3], the introduction of a precise formalism (H [5], the 
reformulation and approximate solution of the combinatorial problems in- 
volved 10 17| . Although the basics are by now quite well understood, there 
still are surprises in store. One example are the structures that were found 
in a computer analysis of the spectrum of states lfl8l[T9| . 

The calculation of the entropy of a non-rotating black hole in loop quantum 
gravity boils down to a rather complicated combinatorial problem. It can 
be treated to a very good approximation in an asymptotic regime [6, 7], and 
proportionality of entropy to area has been established. 

In the present paper we will develop a model that drastically simplifies 
the technical aspects of the entropy calculation while - as our results will 
show - retaining many of the qualitative features of the actual situation. In 
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particular, the combinatorial problem for our model can be solved exactly, 
so any question that one may have about it can be answered with relative 
ease. 

The simplifying assumption that we will make is a rather obvious and sim- 
ple one. One of the hallmarks of loop quantum gravity is a complicated, 
non-equidistant area spectrum. In particular, the area eigenvalues for a 
non-rotating isolated horizon of a black hole are sums of numbers Aj, 

Aj = 87ryl£Vj(j + 1). j G N/2, (1) 

where y is the Barbero-Immirzi parameter, and lp the Planck-length. The Aj 
are obviously not equidistant, however they become approximately equidis- 
tant for large j . Our approximation in this paper consists in using 

Aj = STtylp + (2) 

instead of (QJ, i.e. effectively changing the area operator of the theory. One 
may interpret 10 as the first two terms in the series^] 

v/JITTTT = j + 1 ~ Wttt ~ 1 6(2jV t )3 + • • • < 3 > 

We are certainly not the first use an approximation like this. For example, 
it has been used in [6] to give bounds on the Barbero-Immirzi parameter. 

We should also point out that it has been argued [8] that a very similar 
equidistant spectrum, 

Aj = 87ryl|j j € N/2, 

does arise in loop quantum gravity, upon quantizing the area of a non- 
rotating black hole following an alternative route. We will give results on 
the entropy for this modification of the area spectrum in the appendix. 

The nice thing about the approximation (O is that it simplifies the calcula- 
tion of black hole entropy in the theory tremendously. We can easily cal- 
culate the generating function corresponding to the combinatorial problem 
of enumerating the horizon states. From the generating function, a lot of 
information can then be obtained, as we will demonstrate. 

^ake + 1 } + x, expand around x = 1/4 and evaluate at x = 0. 
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Because of this simplicity our model may be useful, for example to do a first 
quick check on some hypothesis, before attempting to check it for the actual 
system. We demonstrate this by studying - and ruling out - some admit- 
tedly far fetched proposal about describing rotating black holes within this 
formalism (Section |2.3D . 

The paper is organized as follows: In the following section we calculate the 
generating function for the problem and derive various results about the 
asymptotic growth of the number of states and thus the entropy. In Section 
[3] we will show that our results parallel those of [7] for the full spectrum, 
and discuss some ramifications. In an appendix we give results for the 
modified area spectrum A'. 

2 Counting 

In the literature on the subject, slightly different things have been counted 
when calculating black hole entropy in loop quantum gravity (6l [TO]. HH 
[12). This has to do with the fact that one has to distinguish between bulk- 
and boundary-stated and this distinction is not entirely trivial. In practice, 
there arise two different way to count the entropy, and both lead to the 
same results on a qualitative level. 

While we do not want to commit ourselves to either of these ways to count 
on physical grounds here (we may have to say more about this elsewhere), 
we still restrict to only one way to count in the present article. This is merely 
to keep the presentation straightforward. We do not see any problem to 
extend our results to the alternative way of counting C9l [l0l[l2| . 

We will follow the definitions of [5, 6]: Let us call the number of surface 
states with an area smaller or equal a N< ue (A) (the superscript 'true' is 
meant to indicate that this is with respect to the actual area spectrum CD). 
N< ue (A) can be obtained |5j by counting ordered sequences (bOi of integers 
bi modulo k which sum to zero and satisfy certain additional requirements, 
namely: There exist sequences (rrii)i, tru G Z*/2 and ji G N*/2 such 

2 Barring some sort of holography (for which there is currently little evidence in LQG), 
there are infinitely many different bulk states for a black hole of a given area, so counting 
those does not even make mathematical sense. 



3 



thall 

b| = — 2mi mod k, and mt € {— ji> ~h + 1 > • • • > k) 
as well as 

87tyl^ Vii(H+1) < A. (4) 

i 

According to the philosophy laid out in the introduction, we will just change 
the area spectrum, and keep all else unchanged. We will denote by N<(A) 
the number of sequences (b-Oi of integers b^ modulo k which sum to zero 
and such that there are (mJi, and (jOt as above, except for that we ask 



instead of ([4]). In |6| it was shown that the definition of N< ue is equivalent 
to a much simpler one: N< ue is the number of ordered sequences (rat)t, 
mi G Z*/2 such that 

Y_ m t = and 87ryl£ VlmtKlmil + 1) < A. 

i i 

The same arguments can be applied to N< (A). It is easy to see that it is the 
number of ordered sequences (mOi, vrii € Z*/2 such that 



Y_ mi = and Snyl^ Y_ (Wl + \ ) ~ A ' 



(5) 



Let us also define N (A), the number of such sequences that satisfy (0 with 
'<' replaced by '='. 

It was realized in 017) that the counting problem can be simplified by im- 
plementing the two conditions of 10 in separate steps. We will follow this 
strategy and define 

j(rrii,m.2, . ..), tth e Z* : ^~ m; = j,^(|rru| + 1) = a| 



N(a,j) 
Similarly we define 



N<(a,j)=^N(i,j). (6) 



i=1 



3 k = A/47rylp has to be integer - it represents the level of the Chern-Simons theory on 
the horizon (5). 
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Note that N<(A) = N^A/KTrvl 2 ,), 0) etc. 

A useful way to think about the counting problem for N ( a, )) is the follow- 
ing: N (a, j) is the number ways to move, in an arbitrary number of steps, 
on the integer lattice Z, from the point to the point j, such that the total 
length of the path, plus the number of steps, is a. 

The numbers N(a, j) obey a recursion relation similar to the ones given in 
|7| . It is simple to calculate N ( a, j ) for low a using a computer. Here are the 
first few values: 

00000000001 0000000000 
000000000000000000000 
0000000001 01 000000000 
000000001 0001 00000000 
00000001 1 0201 1 0000000 
0000001 20202021 000000 
000001 31 232321 31 00000 

0000 1 43262626234 1 0000 

0001 56396989693651 000 
1 6 10 6 12 14 12 18 12 18 12 14 12 6 1 6 1 

1 7 15 12 1 6 26 20 32 25 34 25 32 20 26 1 6 12 15 7 1 

1 8 21 22 23 42 38 50 53 54 58 54 53 50 38 42 23 22 21 8 1 

where j runs horizontally from -10 to 10 and a vertically from to 11. 
We will now compute the generating function 
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G{g t z)=Y_ L N (a,j)gV. 

a=0 j=— a 

To that end, we refer back to the description of N (a, j) in terms of paths on 
Z. Consider paths with just one step. There is just one such path from to 
j and it has total length |j|. Hence the one-step generating function is 

OO / 1 \ 

G 1 (g,*) = g£(gz)-+(i) = g 2 ( — + j^—) . 



The generating function for paths with n steps is just Gy, and thus we get 
for the generating function for our problem of interest (i.e. paths with arbi- 
trary many steps) 

„ g 2 (z 2 -2gz+l) 

G g,z ) = Y_ Gi g,z n = 7 MU „ \ t±7 J v 

^-j (g + 1) (2zg 2 - (z 2 + z+ 1) g + z) 

Because the N< (a, )) are partial sums of the N ( a, j ), the generating function 
G< (g, z) can be obtained as IIT4| 

oo a -. 

G<(g,z) = yy_ N<(aJ)g Q zi = G(g,z). (7) 

a=0i=-a 3 
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These generating functions contain information about the counting prob- 
lem in a very compact and accessible form. In the following we will extract 
some of this information. 



2.1 The asymptotics of N ( a , ) and N < ( a, ) 

The physical states of the black hole horizon (41 [SJ correspond, in our sim- 
plified model, to the states with ) = 0. Therefore the numbers N(a, 0) and 
N< (a, 0) are of special interest. We will calculate their generating functions 

G'' =0 '(g), G< _0 '(g) and asymptotic behavior. 

The generating function G'' =0 '(g) is the coefficient of z° in G(g, z), 
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G (j=0) (g) = T ^o -G(g,z)dz. 
zm J c z 

where C is a certain contour. Poles of G ( g , z) /z are 

. -2g 2 + g± v / 4g 4 -4g3 + g2-2g + 1 -1 
z = 0, z± = 

2g 

with residues 

g 



Res^Glg.zVz) = - 
Res z± (G(g,z)/z) =± 



g i l 

d-g)g 



(g + 1)\A9-1)(2g-D(2g 2 + g + 1) 
Choosing the contour C around zq and z + gives 

G«=°'(g) = ; (1 ~ 9)9 ^ (8) 

(g + i)^(g-i)(2g-l)(2g 2 + g + l) g + l 

= 2g 4 + 2g 6 + 6g 7 + 8g 8 + 12g 9 + 34g 10 + 58g n + . . . . 

According to 0, the generating function for N<(a, j =0) can be obtained 
as 

Gr°'(g) = T~ — G (j=0) (g) 

i - g 

= 2g 4 + 2g 5 + 4g 6 + 10g 7 + 18g 8 + 30g 9 + 64g 10 + 122g 11 + . . . . 

Let us now look at the asymptotic behavior. The coefficients in the series 
expansion of the second term in ((HJ) are constant, so they do not contribute 
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Figure 1: ln(N (a, 0) ) and ln(N< (a, 0) ) (dots) and the corresponding asymp- 
totic result (solid line) 

at all to asymptotic growth, and we can focus on the first term. Its singu- 
larity at g = 1 /2 is the one closest to the origin and thus we suspect that 
ln(N(a, 0)) ~ ln(2)a to leading order. Since the singularities are algebraic, 
one uses Darboux's lemma (ex. [14]) to verify this. It states that for a function 
of the form f (z) = v(z)(l — z)P with (3 ^ N and v analytic in a region con- 
taining the unit disc, expanding v around 1 gives an asymptotic expansion 
of the function. In particular 

[z n ]f(z) =v(1)[z n Kl -z) p + 0(u- p - 2 ) =v(1)( n "^ -1 ) +0(u-P- 2 ) 

where we have introduced the notation [z n ](. . .) for the coefficient of z n in 
the series expansion around zero. Applied to our case we find 

[g Q ]G«=°'(g) =2 Q [ W Q ]G«=°W2) 

i w-i/1 — w/2 

= 2 Q fw Q lfl -w - s V 



[yv + 2)y/] +z/2 + z 2 /2 



\ 2a fa-y2\_ 1 2 a 



6 \ a J 6v / 7tv / a 
to highest order. An almost identical calculation gives 

[g Q ]G»= 0) (g) 



,u lq (j=o) f < 1 2Q 
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These approximations are compared with the actual N(a,0), N<(a,0) in 
Figured] 

2.2 The asymptotics of unrestricted states 

Our next task is to compute the number of states without taking into ac- 
count the restriction j = 0, 

a a 

T(a)= £ N(aJ), T<(a) = N<(a,j). 

}=-a j=-a 

The generating function for T(a) is simply G(g, 1 ), 

G(g,1)=- - 2 f = 2g 2 + 2g 3 + 6g 4 + 10g 5 + ... 
2g z + g-l 

= l^(2(-l) Q + 2 cl )g a 
a=1 

so that we find T(a) = (2(— 1 ) a + 2 Q )/3 and hence and hence ln(T(a)) ~ 
ln(2)a. For T<(g) we find 

G<(g,1)=-7= rM- T -=2g 2 + 4g 3 + 10g 4 + 20g 5 + ... 

(1 - g)(2g 2 + g - 1) 

5 a=1 

whence T<(g) = (-3 + (-1 ) Q + 2 a+1 )/3 and ln(T<(a)) ~ ln(2)a as well. 

2.3 Asymptotics in both variables 

Finally we take a look at the joint asymptotics for a and ) large and com- 
parable. This is in part motivated by the following curious observation: In 
|7|, the (a, j) asymptotics is calculated to be 

i 2 

lnlN(aJ) ~ cia + C2— 
a 

where Ci and Cz are certain numerical constants of order one. This is rem- 
iniscent of an expansion of the Smarr formula for the area of a Kerr Black 
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hole 




A(MJ) =8ttM z | 1 +A/1-T73 I = 167tM 2 -47t^ + ... 

if one identifies a with M 2 and and ) with J. ) is bounded by a (U so this 
identification would not lead out of the range of allowed spins for a Kerr 
black hole. It should also be noted that a similar suggestion has been before 
[T3|. So, could it be that the states with ) ^ describe rotating black holes? 



There are a lot of reasons to doubt that, including that the meaning of a 
really should be area, not mass squared, that states with j / are un- 
physical according to the framework of |4], and that there is a sophisticated 
treatment of rotating black holes in loop quantum gravity Ifl5| that works 
quite differently. In fact we find that the asymptotics we calculate do not 
match this hypothesis at all, as follows: 

Determining the asymptotics of multivariate sequences using generating 
functions is not a very well known subject and can become quite technical. 
There is however a beautifully developed general theory that one can rely 
on (see for example [16J for the generic case). We will use these results 
in the form presented in [17J. The upshot is that the asymptotics of the 
multivariate sequence in a certain direction is governed by one or more 
critical points, certain singular points z of the generating function. In the 
case of a two variable generating function G(g,z) = I(g, z)/T(g, z), critical 
points are given as solutions z = ( g* , z* ) to 

9 3 
J(g,z)=0, n 2 g — J(g,z) =U!Z— J(g,z) 
og oz 

where (ui , u 2 ) G Z 2 gives the direction in which the asymptotics of the 
sequence is be taken. J and I have to satisfy several properties, for which 
we refer the reader to Jl7|. Which of the critical points actually contribute 
to the asymptotics can be determined by a straightforward but rather te- 
dious analysis which we will circumvent here. If a critical point (g*,z*) 
does contribute to the asymptotics of the sub-sequence c n = c nin , n2n of 
the multivariate sequence c m>n / it does so by a factor g*T nrLl z^ n2TL to high- 
est order. 

We are interested in the asymptotics of N(a, )) for a := j/a constant. With 
the help of mathematica we compute the critical points for the problem at 
hand. There are eight critical points, all depending a. Instead of doing a 
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Figure 2: ln(N(a, j)) (dots) and the asymptotic result (solid line) compared 
for two different values of a 



lengthy analysis as to which of these solutions contributes we simply pick 
the solution with the right limiting values at a = 0, 1 . It reads 

(a - 1 ) ((Y + 1 1 )a 2 + (Y + 4) a + (8a 2 - X(Y - 1 )) a - XY) 
9 * ~ 4a(a+1H-(Y-1)a 2 + X + Y) 

— (Y — 1 )a 2 + X + Y 



2(1 -a 2 ) 
where we have used the abbreviation 



x=v w-ii«» + 4 , y= ,r (5a2 +2X r 31 

V d-« 2 ) 

The asymptotic behavior is thus given by 

ln(N(a,aa)) ~ -ln(gja - ln(zjaa. (9) 

This asymptotic formula works quite well. For comparison we have plot- 
ted (O and the actual values for two different values of a. 

What about the Smarr formula? We have found that, to highest order, 
ln(N(a,j)) is indeed a function of the form af(j/a). However, the func- 
tion f is very complicated and not what one would expect if the suggested 
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Figure 3: Comparison of Ln(N(a, ))) (dots) and the function s from < TTQ]> for 
a = 100 

interpretation were correct. Moreover the numerical fit is quite bad, as Fig- 
ure |3] demonstrates. There we have plotted the numeric result in terms of 
a = j/a (and appropriately normalized), together with the function 

s(a) = 1 + ^J^ - a 2 . (10) 

Thus we have ample reason to throw out the hypothesis that we stated in 
the beginning of this subsection. 

3 Conclusions 

Let us summarize our results and compare them to what is known [7] about 
the counting for black holes using the correct area spectrum To start 
out, all the statements made in |7j for the full spectrum are also borne out 
in our model on a qualitative level. In particular we see: 

• The number of states grows exponentially with area. 

• The number of (un-physical) states with ) arbitrary grows with same 
rate as that of the physical states (j = 0) to highest order. 

• The highest order growth of N< (A) and N (A) is the same. 
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• The next to leading order term in the logarithm of the number of j = 
0-states is —1/2 In (a). 

This confirms the present model as a nice and simplified laboratory' for 
questions about the full theory. Vice versa the present note can be read 
as giving confirmation of the results |7J. We should however note that 
all the numerical factors appearing in our model are different from (al- 
though close to) the ones for the correct spectrum. To give an example, the 
Barbero-Immirzi-parameter obtained in [7 J is 7m ~ 0.24 whereas the one 
that would result if our model were correct is ln(2)/7t ps 0.22. This is not 
so surprising if one remembers a result from [6] that shows that states with 
m higher than 1/2 certainly contribute substantially to the overall count- 
ing, however their numbers are increasingly suppressed with increasing m. 
Thus, among all the states to be counted, the ones for which our approxi- 
mation 0) is relatively bad, are most numerous. 

It is curious that the coefficient of the ln(area)-term in the entropy seems to 
be very robust, as it is -1/2 in the present framework as well as in Q, flHfl, 
and elsewhere. 

There is one last point worth mentioning, one that on first sight seems triv- 
ial: Since the spectrum is equidistant, the number of black hole states, and 
hence the entropy will grow in discrete steps. The height of the steps may 
vary with area, but the size of the steps is always 47rylp. Surprisingly, a 
similar (but certainly more complex) behavior has been observed |T8l [19| 
for the state counting in the full theory. One might at first think that this is 
not an accident. After all, the ladder' observed in IT181IT91 may perhaps be 
understood as a perturbation of the half integer steps in the present model, 
much as the full spectrum ([l} can be viewed as a perturbation 0) around 
the equidistant one. We will investigate this point in much more detail in a 
forthcoming paper (20|. The upshot is that the explanation of the ladder' is 
not that simple. It is related to properties of the area spectrum, in particular 
also to (0, but in a rather involved way. 
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A An alternative modified area spectrum 



The form of the modified area spectrum (0 in the main text was chosen 
so as to approximate the actual spectrum arising in loop quantum gravity 
as good as possible while affording drastic simplification by being equidis- 
tant. There is another equidistantly spaced modification of the area spec- 
trum that has been considered in the literature. It differs form §2$ by re- 
moving the constant 1/2, 

Aj-STtyl^j. (11) 

It can be argued that this spectrum does arise in loop quantum gravity 
upon quantizing the area of a non-rotating black hole following an alterna- 
tive route [8 j@ In this appendix we will derive results for the asymptotic 
behavior of the entropy obtained with the spectrum ((TT)) . We will however 
be more terse than in the main text. We will get rid of all units and some 
constants by considering the quantity 



N(a,j) 



(mi,m 2 ,...), mi£Z\{0} : Y_ m i = i> Y- l m il = a 



A useful way to think about this is the following: N(a, j) is the number of 
ways to move, in an arbitrary number of steps, on the integer lattice Z, from 
the point to the point j, such that the total length of the path is a. 

The numbers N(a, j) obey a recursion relation similar to the ones given in 
|7[. Up to simple changes of variables, they are A035002 in (21|, and they 
are closely related to the combinatorial problems in |22j|l 

It is simple to calculate N(a, j) for low a using a computer. Here are the 



4 This alternative quantization is only possible for the black hole horizon. For other sur- 
faces, only {T) applies. 

5 See especially Section 7 in |22| where the object under study is essentially N (a, 0). 
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where j runs horizontally and a vertically. 
The one step generating function for N is 

G l(g ,z) = £>r+(£r=--V-^. < 12 > 

z — - V z/ q — z qz — I 

n=l 3 ° 

The generating function for paths with n steps is just G^, and thus we get 
for the generating function for N 



nt \ \~ in ( un g(-z 2 + 2gz-l) 

G(g,z)=X(G l(g ,z)) =- 3zg2 _ 2z 2 Q _ 2g i z 



A.l The asymptotics of N ( a, 0) 

The generating function for N(a, 0) is (formally) the coefficient of z° in 

G(g,z), 



G (j=0) (g) = ^o -G(q,z)dz 



2m 



1 



where C is a certain contour. Poles of G ( q , z) /z are 



3g 2 ± V9g4-10g2 + 1 +1 
z± = , z = 

4g 

with residues 

1 1 - g 2 

ReSzJGtg.zVz) = --, Res z± (G(g, z)/z) = ± — = 

2 2 v / 9g 4 - 10g 2 + 1 

Choosing the contour C around zo and z + gives 

i-g 2 1 _ Tg 2 ^ i 



G ( ^ 0) (g) 



2VV"10g 2 + l 2 2V(g-l/3)(g + l/3) 2 
2g 2 + 14g 4 + 106g 6 + 838g 8 + 6802g 10 + O (g 11 ) 
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The singularities g = ±1/3 suggest that ln(N(a, 0)) oc ln(3)a to leading 
order. 'Darboux's lemma' confirms this. We find indeed 

[g Q ]G (i=0) (g) = [v Q/2 ]G (j=0) (v^) = ^9 Q/2 [u Q/2 ](1 -u)- 1/ V9=u 



oc 



W a/2 - 1/2 W4i3< 



a/2 J s/tl sfa. 



to highest order. 



A.2 The asymptotics of unrestricted states 

Our next task is to compute the number of states without taking into ac- 
count the restriction j = 0, 

a 

T(a) = Y_ N ( a 'i)- 

j=-a 

The generating function for T( a) is simply G(g, 1 ), 



_ 2 9 
1 

so that we find T(a) = 2 • 3 a_1 and hence ln(T(a)) oc ln(3)a. 



G(g,D = T -V=2>"3 Q - 1 g Q 



A3 Extremal states 

Out of curiosity, let us also calculate the number of 'extremal configura- 
tions', N(a, a). It this is easiest by going back to ((12)) and modifying it 
appropriately. In this case we only need one variable and we have 

oo 

G 1 (g) = ^g n = T ^ 

whence 

oo n 
n=1 v VJ y a=1 

(This is the well known formula for the number of ordered partitions of an 
integer a.) Hence N(a, a) = 2 a_1 and ln(N(a, a)) oc ln(2)a. 
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Figure 4: ln(N(a, j)) (dots) and the function s from (Tl3] > compared for two 
different values of a 

A.4 Asymptotics in both variables 

Finally we take a look at the asymptotics of N(a, j) for a := j/a constant. 
With the help of Mathematica we compute the critical points for the problem 
at hand. There are four critical points, all depending a. Instead of doing a 
lengthy analysis as to which of these solutions contributes we simply pick 
the solution with the right limiting values at a = 0, 1 . It reads 

_ -L + 3a + 4 VI - a 2 _ /(L — a)a + 2 

9 * ~ 3(a+l) v^2a 2 + 2La + 4' " V 2 -2a 2 

where we have used the abbreviation 

L = ^4 - 3a 2 

The asymptotic behavior is thus given by 

ln(N(a, aa)) oc — ln(g ! „)a — ln(z*)aa. 

What about the Smarr formula in this case? We have found that, to high- 
est order, lu(N(a,j)) is indeed a function of the form af(j/a). However, 
the function f is very complicated and not what one would expect if the 
suggested interpretation were correct. However, the situation is not catas- 
trophic either: In Figure |4] we have plotted the numeric result in terms of 
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a = j/a (and appropriately normalized), together with the function 



s <x 



1 + \/l - a 2 . (13) 



One sees that the numerical correspondence is not bad, but certainly not 
convincing. 
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